Filament depolymerization by motor molecules 
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Motor proteins that specifically interact with the ends of cytoskeletal filaments can induce filament 
depolymerization. A phenomenological description of this process is presented. We show that under 
certain conditions motors dynamically accumulate at the filament ends. We compare simulations 
of two microscopic models to the phenomenological description. The depolymerization rate can 
exhibit maxima and dynamic instabilities as a function of the bulk motor density for processive 
depolymerization. We discuss our results in relation to experimental studies of Kin- 13 familiy 
motor proteins. 
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Many active processes in cells are driven by highly spe- 
cialized motor proteins which interact with filaments of 
the cytoskeleton. Important examples are cell locomo- 
tion, cell division and the transport of organelles inside 
the cell 01 . Cytoskeletal filaments are linear aggregates 
of proteins, for example actin and tubulin. Actin fil- 
aments and microtubules are dynamic and can rapidly 
change their lengths by addition and removal of subunits 
at the ends P, I2 ■ Filaments show a structural asym- 
metry, which provides a direction for motion and force 
generation of bound molecular motors. These proteins 
are able to transduce the chemical energy of a fuel which 
is ATR to mechanical work while interacting with a fila- 
ment 

In addition to generating forces along filaments, mo- 
tors can also interact with filament ends where they may 
infiuence the polymerization rate and thus the filament 
length. Examples are provided by the members of the 
Kin-13 subfamily of kinesin motor proteins A 
particular example is the mitotic centromere-associated 
kinesin (MCAK) which regulates the length of micro- 
tubules during cell division Q . In the course of cell divi- 
sion, the chromosome pairs are separated by the mitotic 
spindle. In this process, shortening microtubules gener- 
ate forces pulling the chromosomes towards the oppos- 
ing poles of the cell. MCAK is localized at the micro- 
tubule ends which interact with chromosomes and it 
has been shown that it induces depolymerization of mi- 
crotubules 0- In vitro assays and single molecule stud- 
ies have shown that MCAK accumulates at both ends 
of stabilized microtubules and induces depolymerization 
at a rate which depends on the bulk motor concentration 
while at the same time MCAK molecules do not generate 
directed motion along microtubules 7]. 

In this paper, we discuss the dynamics of motor 
molecules which induce the shortening of the ends of fila- 
ments to which they bind using both a phenomenological 
description and more microscopic models. For simplicity, 
we consider one filament end and use a semi-infinite ge- 
ometry. The density of bound motors at a distance 2; > 



from the depolymerizing filament end is denoted p{x). 
Here, we use a reference frame in which the depolymeriz- 
ing end is located at a: = for all times. Motors occur in 
the bulk solution at concentration c. They bind to and 
detach from filaments with rates Wac/pmax and cjd, re- 
spectively, where pmax is the maximal density of motors 
for which binding sites on the filament saturate. Bound 
motors diffuse along the filament with a diffusion coeffi- 
cient D and may also exhibit a directed average motion 
with velocity vq ■ Note, that vo in general depends on the 
density of motors p Hfl, El El E3 ■ 
along the filament then obeys 



The density profile 



dtp + dxj = WaC 1 - 



(1) 



The current of motors is given by j = —Ddxp — vp. Here, 
V = Vo + Vd is the total velocity of motors with respect 
to the filament end, with Ud > denoting the depoly- 
merization velocity. It is related to the rate f2 of subunit 
removal from the end by Wd — ^a/N, where a is the size 
of a subunit and N the number of protofilaments in the 
filament. 

We assume that the rate of filament depolymerization 
is regulated by motors bound to the filament end. There- 
fore, the rate of subunit removal is a function Q{po) of 
the motor density po = p{x = 0) at the end. It is useful 
to systematically expand in powers of po 



Vt{po) = + ^ipo + ^2pI + 0{pI) 



(2) 



Here, we have introduced the expansion coefficients fi^. 
The subunit removal rate in the absence of motors Slg 
in general depends on buffer conditions. In situations 
where filaments are stabilized, fJo = 0. For motors which 
induce filament depolymerization, f2i > 0. Since the rate 
f2 saturates for large densities, typically ^2 < 0. 

The description is completed by specifying the bound- 
ary conditions at x = and for x — » 00. At a; = 
the current j{x — 0) at the filament end equals the net 
rate J at which motors attach to the filament end. Since 
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motors attached to the end induce depolymerization, the 
rate J{po) is a function of the motor density po at the 
end and also depends on buffer conditions. Again, we 
express J by an expansion in powers of po- 



JiPo) = Jo + JiPo + J2P0 + C(Po) 



(3) 



Here, Jq is the rate of direct motor attachments to the 
end. The coefficients Ji and J2 characterize how interac- 
tions between motors and the filament end influence the 
detachment rate of motors. If w > 0, motors typically 
detach from the end and thus J < 0. 

Finally, for large x we require that the density p 
approaches the equilibrium value of the attachment- 
detachment dynamics 



Po 



(4) 



If a motor bound to the end removes a filament sub- 
unit, it may fall off the filament with this subunit or 
it may stay bound to the filament. The tendency of a 
motor to stay attached while removing subunits can be 
described by its processivity. In our phenomenological 
description, we define the effective processivity 



PcS = 1 



J{po) - J(0) 



nipo) - n{o) 



(5) 



It differs from the processivity of an isolated motor due 
to collective effects resulting from interactions between 
motors at the end and thus depends on po- If Pcft < 0, 
more motors detach from the end than subunits and the 
motor-induced depolymerization is non-processive. How- 
ever, if < Ppff < 1, a given motor can remove more than 
one subunit 17]. 

For simplicity, we ignore the density dependence of vq- 
In this case, Eq. approaches for large times t, a steady 
state given by p = Poo + (Po ^ Poo) exp[— x/A] The char- 
acteristic length is 



A 



2D 



n 



1/2 



(6) 



The steady state value of po is obtained by inserting this 
expression in Eq. In the following, we consider the 
case where Vd 3> |wo| and the spontaneous velocity vq 
can be neglected. Indeed, experimental observations of 
MCAK show that Vd ^ \vo\ 7]. Depending on param- 
eters, motors either accumulate or deplete at the fila- 
ment end, see Fig.Q] Accumulation at the end occurs for 



(a) 

fl2 > , where 

j-j(a) _ ^1 + J2Pn 
^ Poo 



^0 + JlPn 



JoPn 



PIo 



•(7) 



Motor accumulation can exhibit a reentrant behavior as 
a function of increasing bulk motor concentration, see 
Fig.d 
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FIG. 1: Regimes of accumulation and depletion of motors 
at the shrinking filament end as a function of the coefficient 
= ^2Pm!L^I (D'j^a)^^'^ and the bulk motor concentration 
c* = u^c/uodPi-niL^ for Jo = f2o = J2 = and JiKDuodY^'^ = 
— 0.1, r2i/(_Da;d)^''^ = 0.5. Accumulation occurs above the 
solid line, depletion occurs below. The grey area indicates 
the region of physical interest. Above this area, Pcff (Pmax) > 
which is forbidden by steric exclusion of particles. Below this 
area, filaments polymerize at high motor concentration. 



Our phenomenological description reveals that motors 
can dynamically accumulate at the filament end even if 
their binding affinity to the end is not larger than in the 
bulk, U,a = 0. In this process, motors which bind along 
the filament are subsequently captured by the retracting 
filament end. Dynamical accumulation of motors is a 
collective phenomenon and requires a sufhciently large 
effective processivity. 

In order to obtain a physical picture of the microscopic 
events which influence the effective processivity as a re- 
sult of crowding, we extend discrete stochastic models 
for motor displacements along filaments [lol llll IT^ 
to capture subunit removal at the ends. Motors are rep- 
resented by particles which occupy discrete binding sites 
indexed by i = 1,2,3,... arranged linearly on a fila- 
ment which consists of a single proto-filament (N=l), see 
Fig. 121 Here, i — \ denotes the binding site at the fila- 
ment end. Each site is either empty (n^ — 0) or occupied 
ijii — 1). Particles move stochastically to neighboring 
empty sites with rate uj\-i in both directions. Here, we 
assume again that can be neglected as compared to 
Vd- In addition, particles attach to and detach from the 
lattice with rates WaC and Wd, respectively. The rates of 
particle attachment and detachment at the end site i ~ 1 
differ from the bulk rates and are denoted fiaC and fid. 

Wc first consider the situation where only one particle 
is present. If this particle is not bound at the end, the end 
is stable. If the particle is bound ai i — 1, this subunit 
is removed from the filament with rate fj. This process 
can occur in two different ways: (i) with probability p 
the particle stays bound to the new filament end after 
the first subunit is removed; (ii) with probability 1 — p 
the particle detaches from the filament together with the 
removed subunit. If p is close to one, a single particle is 
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FIG. 2: Discrete model of motor-induced filament depolymer- 
ization. Motors attach to empty sites at rate cDaC and detach 
with rate tJd- The hopping rate to free neighbouring sites is 
Cuh- Occupied sites are removed from the end with rate fi. 
With probability p, the particle remains attached to the end 
when a subunit is removed. 



processive and can repeatedly remove subunits from the 
end without falling off. 

Processive removal of a subunit requires simultaneous 
interaction of a motor with the end and the adjacent 
subunit. Therefore, this process is affected by the pres- 
ence of other particles bound to the filament near the 
end. fn particular, if site i = 2 is occupied while sub- 
unit i = 1 is removed, the new end site is already oc- 
cupied and the particle at i = 1 cannot stay attached. 
We distinguish two cases with different behaviors in this 
situation. Model A describes the case where subunit re- 
moval by a motor requires an empty adjacent binding 
site, while model B corresponds to the case where the 
rate of subunit removal is independent of the occupation 
of the adjacent site. In model A, the probability per unit 
time to remove the end if rii = 1 is fl^ — $7(1 — pn2). 
Here we assume that the processivity characterized by p 
is unaffected by the occupation of the neighboring site. 
Crowding at the end obstructs cutting and reduces the 
rate of subunit removal. In model B, an occupied adja- 
cent site will reduce the processivity of a motor but does 
not affect the depolymerization rate, i.e., Cl^ — ft. 

We can represent the dynamics of the system by a Mas- 
ter equation for the probability P{ni, t} to find a config- 
uration of lattice occupation (ni,n2, ..) at time t. This 
leads to expressions for the rate of change of average oc- 
cupation numbers valid for i >2: 



dt 



i^h {{rii+i) - 2{ni) + + uji,c{l - rii) 



At the filament end, i — 1, 

d{ni) 



Ui)) 



(8) 



dt 



= ujh{n2 - ni) + rJac(l - ni) - il.d{ni) 
- {1 - p)n{ni{l ~ n2)) . 



(9) 



Using a mean-field approximation, replacing two-point 
correlators (riini+i) by (rii) (rti+i), we obtain from 
Eqs. (|H1) and © differential equations and boundary con- 
ditions identical to Eqs. CQ-© with p{x = a{i — 1)) = 
{ni)/a. This procedure leads to explicit expressions 
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FIG. 3: The velocity v*^ = v^/a^uJhUJd)^^^ of depolymerization 
as a function of the bulk motor concentration c* = a)ac/tJd ob- 
tained in simulations of model A for different values of the pro- 
cessivity p = 0, 0.9 and 1 (symbols). For comparison the cor- 
responding solutions of the phenomenological equations are 
displayed (lines). Inset: The accumulation of motors, charac- 
terized by the ratio po/poc of motor density at the end and 
far from the end is shown as a function of c* for the same sit- 
uations. Parameter values are = f^ajWd = fid ~ 0.008 a;h 
and Cl — 4 tJh. 



for the values of the coefficients il; and Ji introduced 
above. For both model A and B we find, using this 
approximation, D = a^oJh, = o-Vl, Jq = f2aC and 
Ji — —a{Cli^c + fid + {1 — P)^)- The nonlinear coef- 
ficients fl2 and J2 are model dependent. In model A, 
0.2 = —a^pfl and J2 — whereas for model B, il2 — 
and J2 = —a^pfl. In both models all higher order coeffi- 
cients fin and Jn vanish. 

Fig. 13 displays the depolymerization velocity Vd ob- 
tained in the mean field theory corresponding to model 
A as a function of the bulk monomer concentration c for 
different values of p. For large c the velocity saturates at 
i'd(pinax)j while it increases linearly for small c. For suf- 
ficiently large p, the velocity Vd exhibits a maximum as a 
function of c. Increasing p further, a dynamic instability 
appears where two stable states with different Vd coexist 
within a range of c values. A third unstable state is indi- 
cated by a broken line. Results of stochastic simulations 
of model A are shown for comparison. Mean field the- 
ory and stochastic simulation agree quantitatively except 
in the vicinity of the dynamic instability. Fluctuations 
conceal the dynamic instability present in mean field the- 
ory. The inset to Fig. O shows the relative accumulation 
Pa/ Poo of motors. For sufficiently large p, motors accu- 
mulate as c is increased. We note that in model B, no 
dynamic accumulation of motors occurs and pQ < poo- 

In summary, we have shown that a positive effective 
processivity peff of subunit removal is essential to achieve 
dynamic accumulation of motors at the filament end. 
This effective processivity is a collective effect and re- 
sults from steric exclusion of motors bound near the end. 
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The phenomenological description given by Eqs. Q-® is 
general and valid irrespective of details of the mechanism 
of motor induced subunit removal and of the structure of 
the depolymerizing filament end. We have restricted our- 
selves to effects corresponding to the lowest order terms 
in the expansions of Eqs. (|2Jl and While higher 

order terms could lead to additional effects, our simula- 
tions of microscopic models indicate that these terms are 
unimportant (see Fig. 3). We have focussed on excluded 
volume effects at the filament end and have considered 
the case ?;o = where some of these effects in the bulk 
disappear. For vq ^ 0, the interplay between bulk and 
end excluded volume effects could lead to new phenom- 
ena which will be subject of future work. 

The stochastic models A and B provide a physical 
picture of the cooperativity and processivity of motors 
bound at the end of a single protofilament. Interac- 
tions between motors lead to different rates of subunit 
removal in the two models. In our stochastic simulations 
for iV = 1, the dynamic instability of steady states which 
is found in the mean field analysis is concealed by fluc- 
tuations. We expect that for larger numbers of protofil- 
aments this effect of fluctuations is reduced. Therefore 
a signature of a dynamic instability could reappear for 
microtubule depolymerization leading to bistability and 
switch like changes of depolymerization velocities. In the 
mitotic spindle this instability could be relevant for chro- 
mosome oscillations, which have been observed [l^ . 

Accumulation of motors at the filament end described 
by Eq. (9) can occur as a result of three different mech- 
anisms. Motors can accumulate by directly binding to 
the filament end if they have a higher affinity to the 
end than to subunits along the filament. This effect 
dominates if the total velocity v = t;o + '^'d is small, 
^ 4:D{ujac/ pmax + ^^d)- In this case, the localiza- 
tion length is given by the diffusion length during the at- 
tachment time A ~ -D^/^(wac/pmax + Wd)~^^^. A second 
mechanism of accumulation is given by transport of mo- 
tors to the end with velocity vq » Wd- In the third case, 
motors that bind along the filament are captured by the 
shortening end. This dynamic accumulation mechanism 
dominates for vq ^ and v"^ ^ 4Z?(a;ac/pmax + t^d)- 
The localization length is A D/vd- The first and last 
cases can lead to accumulation at both ends of a filament 
of finite length, while in the second case motors accumu- 
late at one end only. 

Our results can be related to experiments on members 
of the Kin-13 family of kinesins. The depolymerization 
velocity Vd as a function of bulk motor concentration has 
been measured for MCAK and it has been shown that 
MCAK accumulates at both ends 0]- The observed ve- 
locity Vd is consistent with both models A and B since 
it does not exclude the possibility of a maximal velocity 
for intermediate motor concentrations. Accumulation at 
the end suggests that a mechanism similar to model A is 
more likely to be at work. Indeed, experiments indicate 



that a collection of MCAK motors processively depoly- 
merize microtubules consistent with model A. Present 
data cannot rule out a mechanism akin to model B where 
motor accumulation is still possible if the affinity of mo- 
tors to the filament end is high. In future experiments, 
model B could be ruled out if a a maximum of the de- 
polymerization velocity at intermediate motor concentra- 
tion would be observed as suggested by our theory. The 
members XKCMl and XKIF2 of the Kin-13 family, can 
depolymerize microtubules with or without accumulation 
of motors at the end, depending on the conditions under 
which microtubules have been stabilized [^. Further- 
more, it has been suggested that processivity is reduced 
under conditions where motors do not accumulate 0. 
Thus, stabilization of microtubules could influence the 
microscopic mechanisms of collective subunit removal by 
a change in the microtubule lattice structure, leading to 
reduced processivity p or a mechanism similar to model 
B. 

The theory developed here is not restricted to motors 
of the Kin-13 family which interact with microtubules but 
applies in general to associated proteins which regulate 
the dynamics of filament ends. Actin depolymerization 
by ADF/cofilin as well as the polymerization of actin by 
formin are further examples of such processes 0, Il6l| . 
In addition to the conventional action of motor proteins, 
filament polymerization and depolymerization by proces- 
sively acting end-binding proteins are expected to play a 
key role in cytoskeletal dynamics and self-organization. 
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